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We extend the Keldysh formulation to quantum spin systems and derive exact equations of motion. This 
allows us to explore the dynamics of single spins and of ferromagnets when these are inserted between super- 
conducting leads. Several new effects are reported. Chief amongst these are nutations of single S=l/2 spins 
in Josephson junctions. These nutations are triggered by the superconducting pairing correlations in the leads. 
Similarly, we find that on rather universal grounds, magnets display unconventional spin wave dynamics when 
placed in Josephson junctions. These lead to modifications in the tunneling current. 
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PACS numbers: 71.27.+a, 71.28.+d, 77.80.-e 



I. INTRODUCTION 

There is a growing interest in a number of techniques that 
allow detection and manipulation of a single spin. A par- 
tial list includes optical detection of electron spin resonance 
(ESR) in a single molecule tunneling through a quan- 
tum dot and, more recently, the ESR-scanning tunneling 
microscopy (ESR-STM) technique H, 0]. Interest in ESR- 
STM lies in thepotential of detection and manipulation of a 
single spin jsl 0]- an ability which is crucial to spintronics 
and quantum information processing. Much work also ad- 
dressed coupling, feedback effects, and decoherence in a cou- 
pled electronic-vibrational systems, such as nanomechanical 
oscillators and local vibrational modes f^. In particular, spin- 
tronic and quantum computing applications greatly intensified 
interest in Josephson junctions. In a previous publication (Ml, 
four of us studied the effect of the supercurrent on a macro- 
scopic spin cluster (of spin S ^ 1) precessing in the pres- 
ence of a magnetic field when placed in a Josephson junction 
to find new spin dynamics. In j^, these systems were ex- 
amined anew wherein variations in the DC current were pre- 
dicted for all systems harboring a spin of any finite size S. 
In the current article, we complement 0] by studying, for the 
first time, the explicit dynamics of single quantum 5 = 1/2 
spins in Josephson junctions to find new intriguing dynami- 
cal effects for which we provide quantitative expressions. The 
single spin {S = 1/2) dynamics which we study here dif- 
fers significantly from the the large magnetic cluster {S ^ 1) 
dynamics studied in In the current article, we further ex- 
amine spin wave dynamics in ferromagnets when placed in 
Josephson j unctions . 

The analysis of spins embedded in Josephson junctions has 
a long and rich history. Early on, Kulik Il0i1 argued that spin 
flip processes in tunnel barriers reduce the critical Josephson 
current as compared to the Ambegaokar-Baratoff limit ( Till . 
More than a decade later, Bulaevskii et al. Iil2i1 conjectured 
that TT-junctions may form if spin flip processes dominate. The 
competition between the Kondo effect and the superconduc- 
tivity was elucidated in 1 1 3 ] . A nice review of experimental 
works on certain aspects of magnetic nanoparticles in Joseph- 
son junctions is found in 11411 . Transport properties formed the 



central core of many pioneering works, while spin dynamics 
was relegated a relatively trivial secondary role. In the current 
article, we report on exact new non-stationary spin dynamics 
and illustrate how a quantum 5 = 1/2 spin is affected by the 
Josephson current. As a consequence of the Josephson cur- 
rent, spins exhibit non-planar precessions while subject to the 
external magnetic field. As well known, a single quantum spin 
in a magnetic field exhibits circular Larmor precession about 
the direction of the field. As we report here, when the spin is 
further embedded between two superconducting leads, quan- 
tum pairing correlations lead to new out-of -plane longitudinal 
motion, much alike that displayed by a classical mechanical 
top will arise. We term this effect the Josephson nutation. 
Similar effects occur when a ferromagnetic slab is placed be- 
tween two superconducting leads. We outline how transport 
is, in turn, modulated by this rather unusual spin dynamics. 
The coupling of the spin to the supercurrent leads to an effec- 
tive non-local in time interaction of the single spin with itself. 
Keldysh contour calculations illustrate that a non-local in time 
single fermion action is also found in situations wherein the 
single spin is replaced by an Anderson impurity JTsll . As well 
known, in the limit of small hopping amplitudes to and from 
an Anderson impurity, the impurity attains a Kondo like char- 
acter much like that of the single spin which is the focus of our 
attention. Here we consider the origin of this rather generic 
non-locality in time present in the dynamics of a Josephson 
junction. En route to deriving this new spin dynamics we il- 
lustrate that even in the presence of non-local in time inter- 
actions, certain variants of the classical equations of motion 
become trivially exact by virtue of compactness of the spin 
variables. An elaborate extension of these ideas will be de- 
tailed elsewhere lllf 



II. OUTLINE OF THE ARTICLE 

The main goal of the current publication is to report on the 
spin and spin wave dynamics (of single spins and of magnetic 
systems, respectively) in Josephson junctions. 

To achieve this aim, we will initially (in Sections( llIIIIV» 
extend the non-equilibrium Keldysh formalism to address 
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these problems. In Sectionl HIH . we illustrate that even in the 
presence of effective non-local of time interactions of a spin 
with itself (such as those borne by the interaction of a single 
spin with a Josephson current), the equations of motion un- 
dergo a trivial modification. In Section dTvl i. we rewrite these 
equations within the standard Keldysh basis best suited for 
non-equilibrium problems. Sections dilll IIV> may be seen as 
independent extensions of basic facets of the non-equilibrium 
Keldysh formalism for a spin system. 

In Section(|V|i, we apply the rather general formalism de- 
veloped in Sections dIIIIIVt to the specific problem of a sin- 
gle 5=1/2 spin in a Josephson junction (with a time inde- 
pendent potential difference between the two superconduct- 
ing leads). We start, in subsection dV A> . in writing down 
the relevant Hamiltonian of such a Josephson junction har- 
boring a single spin. In subsection JV B> . we briefly high- 
light the natural time scales in the problem- which will in- 
deed come to the fore in the detailed solution which we will 
later expose. In the all-important subsection llV C> . we high- 
light the origin of the effective non-local in time interactions 
of the spin with itself. Here, we integrate out the lead elec- 
trons to find the effective spin only action harboring such non- 
local in time interactions. These non-trivial interactions are 
the reason that we needed to develop and extend (Sections 
(IIIIIIV» the Keldysh formalism to a very general spin sys- 
tem with such interactions. In subsection dV D> . we invoke 
the results of Sections dlllllVl i to the resultant effective spin- 
only action of Section dV CI to write down the equations of 
motion for the spin. In subsectionl lV E> . we solve these equa- 
tions of motion to lowest order in the spin-dependent tunnel- 
ing amplitude. Detailed technical aspects of the solution on 
which subsectionl lV E> dutifully relies on have been relegated 
to appendices B and C. The perturbative solution to the equa- 
tions of motion- the final equations of subsectionlj^^- form 
one of the main core results of the current publication. In 
subsection llV F> . we examine the physical meaning of this so- 
lution of the single spin problem to unearth several new pre- 
dictions for this 5 = 1/2 system. In this subsection, we aim 
to further arm the reader with an intuitive understanding for 
the physical origin of these new effects. Some of these pre- 
dicted effects (and our prediction of nutation in particular) are 
highlighted in Fig.l|3j. In subsection llVH . we examine the be- 
havior of the system for a single spin of magnitude S > 1/2. 
In the large 5 S> 1 limit, we recover our very different semi- 
classical spin (S oo) results of H. 

Next, in Section llVIH . we discuss a variation of the single 
spin problem wherein an AC voltage bias is applied across 
the Josephson junction. Our main result are the predictions 
of specific time dependent spin dynamics displaying an infi- 
nite number of harmonics and new DC lock-in effects. The 
predicted supercurrent in this system is also discussed. 

In Sectionl fVIIH . we examine the problem of a ferromagnet 
in a Josephson junction. In the spin-wave approximation, we 
find that each spin-wave mode displays some of the unusual 
effects predicted in subsectionl lV El IV F> for the single spin 
problem. The predicted spin wave dynamics and associated 
transport (current), are furnished. In Sectionl lDO i. we discuss 
simple extensions of our results to other systems generated by 
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FIG. 1 : The standard Keldysh contour. The times T and T' are taken 
to be — oo and oo respectively. The form of this contour will be heav- 
ily employed in our work when time ordering various spin products. 



a trivial change of geometry wherein at least one of the su- 
perconductors forming the Josephson junction is replaced by 
a planar superconductor In SectionljXj, we write down the 
5^1 equations of motion for general magnets and antifer- 
romagnetic chains. The non-uniform temporal evolution of 
each of the spin-waves is highlighted in the resultant solution. 
We conclude the main text, in Sectionl lXll . by highlighting 
our conclusions. 

In Appendix A, we briefly discuss several experimental 
manifestations of our effect and highlight a proposed exper- 
iment which may verify our predictions. 



III. EXACT SPIN- 1/2 EQUATIONS OF MOTION ON 
KELDYSH CONTOURS 



We start by deriving the equations of motion for a very 
general spin- 1/2 system having two (or more) local and non- 
local spin-spin interactions at different times. In this work we 
employ the non-equilibrium Keldysh technique. Within this 
framework, the spin operators on both up and down portions 
of the (Keldysh) contour of Fig. 1 are normalized and satisfy 
[Su{t), Sd{t')] — 0. In what briefly follows, we will dispense 
with operator formulations and employ a path integral repre- 
sentation. Towards this end, our working horses will be the 
CPi spin coherent variables (z) 

(13 El 

wherein the spins are 

represented by 



S = Sz*CTabZb 



(1) 



(with S the spin magnitude). Here and throughout, we set h — 
1. In Eq.Q, a, h =1, | and we assume an implied summation 
over repeated indices. The vectors Sab are the ah components 
of the three Pauli matrices. The components Za=-\.i code for 
a two component complex spinor subject to the normalization 
constraint, + = 1. By glancing at Eq.Q, we note 
that a knowledge of S specifies the two component spinor z 
only up to a global multiplicative phase. 

As well appreciated, in a spin coherent basis, the Berry 
phases associated with the spin coherent states are the net area 
of the spherical triangle spanned by the spin as it moves on the 
Bloch sphere. The latter may be expressed in the CPi basis as 
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Sserry = i J dt J2a ^l^tZa fivlfl^ . For the benefit of read- 
ers unfamiUar with this formalism, we provide in 1193 a quick 
derivation for this form of the Berry phase. 

We now assume the action contains the single spin term 
~Va J dtSa ■ h describing a single spin in an external magnetic 
field set by h. The parity rja = ±1 is fixed by the direction of 
the contour- ijup = 1, Vdown = We further include a non- 
local in time spin interaction i]aVb J dt J dt' Kab{t, t')Sa{t) ■ 
Sb{t'). The kernels Kab encapsulate non-local temporal de- 
pendence. The generalization to higher order terms is straight- 
forward and leads to no qualitative change. With the Berry 
phase included, the general action 



any 6x{t) leads to no change to the value of Z- the range 
of integration in Z = J Dx e^^ is unchanged. This, in 
turn, mandates that (^)s = fl6ll . Next, we consider 
A'' = z*a^j^ and explicitly illustrate that its expectation 
value vanishes, (^') = 0. To this end, we write the spinors, 
longhand, in terms of real and imaginary components, z* = 
(^fle~*^/Tri ^Re~''-'^im) ^nd the msasuie DzD z* S{\z\'^ —1) = 
DzkDzj^DzleDzlJiKJ^ + KJ^ 



'/ml 



1). Here and in what briefly follows we suppress a uniform 
Keldysh contour index. The expectation value (A*) for each 
value of the spin index, i = x,y, z, is an integral overbilinears 
in z and hence amounts to a sum of integrals of the type 



S = 2iSrja / dtz*dtZa - S dtrjaha ■ z*aza 



/" 1 1 2 2 

lap = / Dz^^Dzjj-^Dzi^^Dzj^ 



+ S^TjaVb dt dt'Kabit,t')z*{t)aZa{t)-z;{t')dzbit').(2) ^(l^^el +\zim\ +\zRe\ " l)^^ 



6S 

Szf3 



iS 



(7) 



Varying the action, 

6S 



Srja[2idtZay{t) - h ■ a^^Zay{t) 



Sz*aAt) 

-T]b [ dt'Kab{t,t')zl^{t')aySZbs{t') -ff^jZajit)^ 

EE Srta (2ldtZa,{t) - H (t) ■ . (3) 



Here, 



Kabit,t') = Kabit,t') + Kbait' ,t) 

Hit) = h + Srja [ dt'Kabit, t')zl^{t')a^szbs{t'). (4) 



Next, we briefly generalize Ehrenfest's theorem to situa- 
tions such as the one of relevance here where a non-local in 
time action is present. A full discussion of this theorem for 
general systems will be presented elsewhere In what fol- 
lows, the expectation value of any quantity A evaluated with 
the action S is denoted by 



DzDz*5{\z\^ - l)Ae 



iS 



(5) 



with Z = / DzDz*5{\z\^ — l)e''^ the associated partition 
function. Similar definitions apply, with a trivial replacement 
of the measure when the action is a functional of one of more 
real fields {xa (i)}- the current context, .Tq, code for the real 
or imaginary parts of the complex spinor components z. Next, 
we note that for any cyclic coordinate x, the expectation value 
of the variational derivative. 



(^)5-^[e -0. 



(6) 



In the above, by the compactness of x, in integrating all possi- 
ble trajectories x{t), the initial and final trajectories are equal 
Xi{t) = Xf{t). This in turn lead to the vanishing expecta- 
tion value given in Eq.(|6} for all non-singular actions. Anal- 
ogously, this result follows by noting that for compact co- 
ordinates, the transformation [x{t) — * x{t) + 6x{t)], with 



Here, the indices a and [3 span the four possible values 
(1 Re, 1 /to, 2 Re, 2 Im). An immediate consequence of the 
vanishing of the expectation value (||) for any cyclic coor- 
dinate X is that all integrals of the form la^fj vanish. An 
inspection of (A^) reveals that the contributions of all inte- 
grals of the type laa cancel identically when i ~ z (the only 
place where integrals of the type Ia=p appear). Similarly, 
{A''') = is^cF'z) = 0. The vanishing {A') = {A'^) = 
imply that their difference. 



0=([ 



SS 



5Zaa{t) 



(7^ z"-'^ — Z* 



5S 



(8) 



where the Keldysh contour index (a) is reinstated. 

Next, we explicitly insert Eq.(|3j into Eq.(|8j. As a conse- 
quence of the SU(2) algebra of the Pauli matrices, we find 
that for each Keldysh contour index a = top/bottom. 



, dSa , 



{H X Sa)s- 



(9) 



Eq.(|5J is none other than the equation of motion for preces- 
sion of the spin S in the instantaneous field given by H of 
Eq.@. We find that such classical equations of motion for a 
nonlocal in time action are exact in the quantum arena. [For 
affectionados of parafermion methods, we briefly note as an 
aside that although throughout we employed the bosonic spin 
coherent path integral representation, a similar result follows 
if the spinors z were Grassmann variables (a net even num- 
ber of permutations of the spinor coordinates are involved in 
proving Eq.(|9j).] The bulk of the paper will be devoted to 
a solution of Eq.Q for different realizations of a Josephson 
junction system. 

We will momentarily dispense with the Keldysh contour in- 
dices. Due to the commutation relations Sx S = iS, although 
the field H contains a piece which is linear in S, the planar 
components of Eq.(|9}, may be reduced for certain problems 
to a linear equation in planar spin components (Si) (i = x, y) 
which then must have the solution 



(5,(0) =C/.,(t)(^,(0)). 



(10) 



3 



LA-UR-04-6501 



February 2, 2008 



submitted to Physical Review B 



We now invoke symmetry constraints. An external magnetic 
field h in the action (Eq.(|2}) lifts the SU{2) spin rotational 
symmetry of the free spin leading in turn to a lower U{1) sym- 
metry of rotations about the external magnetic field axis. Such 
a symmetry is trivially encapsulated by the operator (6) ro- 
tating (5) by an angle 6 about the z (or magnetic field) axis. 
As a consequence, the evolution operator U{t) of Eq.( I10> 
must commute with R^{6). This, in turn, dictates that if the 
solution is in the form of Ea.( llO> . then the time evolution op- 
erator U{t) must have the form 



p{t) q{t)\ 

U{t) - I -q{t) p{t) . 



(11) 



Similarly, due the azimuthal rotational symmetry encapsu- 
lated by R^{0), the expectation value {Sz{t)) must be inde- 
pendent of {Sx{0)) and (5^,(0)). This form will indeed be 
borne out for our full Keldysh problem. 



IV. THE KELDYSH BASIS EQUATIONS OF MOTION 

Within the non-equilibrium Keldysh formalism it is often 
advantageous to apply a simple linear transformation from the 
basis of up and down contour fields to the symmetric and an- 
tisymmetric linear combination of these fields. E.g., for the 
spin 

1 



(12) 



The utility of this basis has its roots in the natural form for 
the various correlation functions- all simply related to the ad- 
vanced, retarded, and "Keldysh" correlators. The subscripts 
"cl" and "qu" of Ea.( ll2> coding for "classical" and "quan- 
tum" suggest an intimate relation to classical and quantum 
Langevin like dynamics. We refer the uninitiated reader to 
excellent texts such as ll20ll . ll2lll where the origin of this link 
is explored in depth. In Ea.( ll2t we trivially generalize this 
change of basis to quantum spin systems. In this basis, when 
taken as operators in Ea. (ll2> [prior to a passage to a path inte- 
gral representation], the spins no longer obey canonical com- 
mutation relations the spins no longer obey canonical com- 
mutations relations (e.g., [Sqm Sd] ^ 0) and are no longer 
normalized {Sup ± Sdown may correspond to a spin-triplet, 
S = 1, or to a spin singlet, S = 0). Thus, we may not directly 
employ the CPi representation in this basis. For the current 
purposes, the equations of motion in this basis may be derived 
from Eq.(|9} for the up and down contour spins. 



= 



-Ka 



4 



^^'^'^^ US 



(13) 




FIG. 2: Magnetic spin coupled to two superconducting leads. 



and 



K., 



-Kd 



Kuu-K^d+Kdu-Kdd 



(14) 



An average over exp[iS'] is implicit in ( )5. As emphasized 
earlier, these are not merely saddle point equations but are 
rather exact. In the above, although the time arguments were 
not explicitly written down, Kap serves as a shorthand for 



V. SINGLE SPIN DYNAMICS IN A JOSEPHSON 
JUNCTION 

A. The system 

Our system is sketched in Fig.|2] It consists of two identi- 
cal ideal superconducting leads coupled each to a single spin; 
the entire system is further subjected to a weak external mag- 
netic field. In Fig.(|2j, ijll.r, denote the chemical potentials of 
the left and right leads, S is a weak external magnetic field 
along the z-axis, and S = {Sx, Sy,Sz) is the operator of the 
localized spin. The wave-functions of our system are super- 
positions of the direct product of states of the left contact, the 
impurity spin, and the right contact. 



(15) 



A tunneling matrix couples these different states. The Hamil- 
tonian of this system reads 



Ti. = Ti.0 + Ti-T, 



l-^ Rka 



acx' • S 



-h.cil6) 



Here, Hl and Hr are the Hamiltonians in the left 
and right superconducting leads, while c]f^^ (cika) cre- 
ates (annihilates) an electron in the lead a in the state 
k with spin a in the right/left lead for i = L/R 
respectively. nL{R) = Efc(p);<T ^Mp)4(p),<rCfc(p),'- + 
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, where we of freedom and Tk denotes time ordering along the Keldysh 



denote the electron creation (annihilation) operators in the left 
(L) lead by c^^ (cka) while those in the right (R) lead by cj,^ 
(Cpcr)- The quantities k (p) are momenta, a the spin index, 
while £k{p},<T ™d Actct' {k{p)) are, respectively, the single par- 
ticle energies of conduction electrons, and the pair potential in 
the leads. In EQ.( ll6t . the components (laa' are entries of the 
three Pauli matrices {(J^a'^ '^aa'^ '^aa')- '^^e current pub- 
lication we consider s-wave symmetry pairing in the super- 
conducting leads. Here, ji is the magnetic moment of the spin. 
With the spin embedded in the tunneling barrier, the conduc- 
tion electron tunneling matrix becomes, spin-dependent |0] 
T = [To 1 + Ti S* • <Tc] . Here To is a spin-independent tunnel- 
ing matrix element and Ti is a spin-dependent matrix element 
originating from the direct exchange coupUng J of the con- 
duction electron spin CTc to the localized spin S. Henceforth, 
we will omit the c subscript. We take both tunneling matrix 
elements (Tq and Ti) to be momentum independent. This is 
not a crucial assumption and is merely introduced to simplify 
notations. Typically, from the expansion of the work function 
for tunneling, 21 ^ J/U, where U is the height of a spin- 



independent tunneling barrier 12211 . A weak external magnetic 
field Bz ~ 100 Gauss will not influence the superconductors 
and we may ignore its effect on the leads. The operator e*"^/^ 
is the (single electron) number operator When the junction is 
linked to an external environment, the coupling between the 
junction and the environment induces fluctuations of the su- 
perconducting phase difference across the junction {4>{t) ). 



B. Physical Time Scales 

The Josephson junction with the spin has two time scales: 
(i) The Larmor precession frequency of the spin tj^ = 
gfisB = h, where g, fiB are the gyromagnetic ratio and Bohr 
magneton of the conduction electron, respectively, (ii) The 
frequency ujj = 2eV, with e the electronic charge, charac- 
terizes the Josephson effect when an external voltage V is ap- 
plied across the junction. 



C. The Effective Action 

Josephson junctions are necessarily embedded into external 
electrical circuits. This mandates that the dynamics explicitly 
depends on the superconducting phase difference (j){t) across 
the junction. The evolution operator is given by the real-time 
path integral 

Z = J D(j)DS cxp[iS] . (17) 



The net action of Eq.l^J is given by S* — [5c 



5spin(5') +5tunnoi(</>, 5)]. The effective action 5tunnei contri- 
bution describes the junction itself. If all external fields are the 
same on both forward and backward branches of the Keldysh 

1, 



contour The label J)^ denotes integration along the Keldysh 
contour as shown in Fig.([0. We first take a partial trace in 
Z over the lead fermions (the bath) to obtain an effective 
spin action. The Josephson contribution to the resulting spin 
action reads 



much in the spirit of Refs. l23l l24l 1251 . For brevity, we set 



The tunnehng Hamiltonian of a phase 



(voltage) biased junction 

T^T = [ToSaa' + TiS ■ CFcrcr'] 

X {A,,, exp(i(/./2) + Al^, exp(-^0/2)) . (18) 

In the presence of a dc voltage bias, cf> = 2eVt. If <f) is treated 
classically (i.e. (p is the same on the upper and the lower 
branches of the Keldysh contour), the contribution cx Tg to 
SS vanishes. The mixed contribution cx TqTi vanishes due to 
the singlet spin structure of the s-wave superconductor The 
only surviving contribution reads 

- ^ £ rft £ dt' [S(i) • aa^fj] [S(t') • as^] 

X (^{TKAc.p{t)As^{t'))e'^^^^ + (A, ^ At, -0))(19) 

where we keep only the Josephson (off-diagonal) terms. The 
spin structure simplifies for the s-wave case: 

Tl i dti dt' [S(i) • S(i')] [iD{t, t')] , (20) 

JK J K 

where the kernel iD{t, t') is dictated by 

{TKA^^{t)A^^{t'))e'^^'^ + {A,^ ^ At,-0). The 
operators A are bilinears in Fermi operators and thus the 
correlator (T/^yl|^(t)A||(t')) will amount to a sum of a 
product of two terms: a product of two normal Green's 
function G and a product of two pair correlators F. Thus, 
generalizing the known effective tunneling action for a 
spin-less junction i23l l24l l25ll to the new spin-dependent 
arena, we obtain 



5tun„ei = -2<t dti dt' a{t, t') + TlS{t) ■ S{t') 
Ik JK 

-2(f) dti dt'P{t,t') -TlS{t)- S{t') 
Jk Jk 

- 4>{t) + 4>{t') 



(71) 



contour (A') then Z ~ TrT^ cxp[- 



TK 



dtHrit)] 



where the trace is over both the electron and the spin degrees 



where ia{t,t') = G{t,t')G{t' ,t) and i(3{t,t') = 
F{t, t')F'< {t, t'). Here, the Green functions 

G(t,<')^-*5](rKCk,(t)4,(t')), (22) 

k 

F{t, t') = -iY^ {Tkc^^ (i)c-ki it')) , (23) 

k 

F\t,t') ^ -^Y.{TKcl{t)cl^^it')) . (24) 
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We now express the spin action on Keldysh contour in the 
basis of coherent states 



^spin = ^ (p dt h ■ S + S' 



WZNW 



(25) 



K 



The second, Wess-Zumino-Novikov-Witten (WZNW), term 
in Ea.( l25t depicts the Berry phase accumulated by the spin 
which we discussed earlier in the coherent spin representa- 
tion wherein it amounts to a kinetic bilinear- the first term of 
Eq Q). In the calculations that follow we replace the spin 
measure DS by the coherent spin state measure DzDz* and 
rely on our derived exact equations of motion. We now per- 
form the Keldysh rotation of Eq. J12t . defining the values of 
the spin and the phase variables. For the superconducting 
phase, we introduce (with notations following Refs. jS^IS) 



(26) 



Within the Keldysh framework, the Josephson current is given 
by 



(27) 



with $0 the unit fluxon (with full units restored, $o — hc/e 
with c the speed of light). With these definitions in hand, the 
tunneling part of the action reads 



'tunnel 



— Sa+ Sr 



(28) 



the pieces together, we find that 



A2 



P^it^t') = -e{t-t') J2 sin m + Ep){t ~ t')] . 

k.p " ^ 

(30) 

The kernel f3^{t — t') decays on (short) time scales of order 
0{h/A). Similarly, 



k,p 



Henceforth, we will often employ the shorthand 
(3^/^[t,t') = P"/^{t - t'). Looking at Eq.lEll, we 
see that the Fourier transform (3^ (cj) vanishes for frequencies 
w < A. This is not so for the retarded correlator 
due to the presence of the theta function. For now, we 
ignore the fluctuations in the superconducting phase and set 
(l>up{t) = 4'down{t) = 4>{i) = ^.Jt with ujj = 2eV (and thus 
X — 0). In this, "classical", limit 

Stunnel^^ j dt j dt' ^"^'{1 , t' S qu{t) ■ S ,l{t' ) ] {t , t' ) 

+ [ dt [ dt'Tfp'' {t,t')S,.^{t) ■ S,u{t')j{t,t'){32) 



withjft,<0 = cos ^'*'t^'*'^ - 



where the normal (quasi-particle) tunneling part is ex- 
pressed via the Green functions = 0{t — t'){a^—a'^) and 



ico) 



where m>(t,i') = G> (t,t')G< {f ,t) 



andia<{t,t') = G<{t,t')G> {f ,t). Similarly the Josephson- 
tunneling part 5^3 is expressed via the off-diagonal Green's 
functions = 9{t - t'){P> - /3<) and /3^(w) = /3> + /3<, 
where iP>{t,t') = F> {t,t')F^> {t,t') and iP<{t,t') = 
F<{t,t')F^<{t,t'). The pair correlators i^<(t,t') are derived 
from F^{t, t') by the interchange of t with t' . In the current 
article, we focus on the interaction between the supercurrent 
and the spin. 

In Eq. J21> . the normal-tunneling part Sa is obtained from 
Sp by the following substitution: f3^/^{t, t') a^/^{t, t'), 
(t){t') -<l)it'), and x(t') -xit')- The Keldysh terms 
(those including f3^ and a^), which normally give rise to ran- 
dom Langevin terms (see, e.g., Ref. L25>] ) are, in our case, sup- 
pressed at temperatures much lower than the superconducting 
gap (T ^ A), due to the exponential suppression of the cor- 
relators (}^{uj) and a^{uj) at < A. 

To obtain we start from the Gor'kov Green functions 



A 



F>\t,t')^iy—. 



,-iEk(t~t') 



(29) 



where the quasi-particle energy Ek = i c^,, 

free-conduction-electron dispersion in the leads. Putting all of 



e|, with Efe the 



D. The Equations of Motion 

With the action at our disposal, we now write down the ex- 
act equations of motions and give a solution, exact to order 
0{T'i). Extracting, in the up-down contour basis, the co- 
efficients, 'Kab{t,t') of the Sa{t) ■ Sb{t') terms in EQ.(l32t. 
constructing Kab{t, t') from Eq.0, and invoking Eq.(I13>. we 
find 



= {-Sa + hxS,i 
dt 

+4|ri|2 / dt'j{t,t')/3''{t,t')Sa{t') X Sa{t) 



+2\Tif / dt'j{t,t')p''{t,t')Sgu{t') X Sa{t) 



+ \Ti\^ J dt'3{t,t')(3''{t,t')Squ{t') X Squ{t))s 

= (^"^c; + X S'c/ + Icl-d + Iqu-cl + Iqu-qu)s- (33) 

The final subscript S serves to remind us that this is the path 
integral average computed with the action S. The various 
subscripts of the integrals / denote the terms that they orig- 
inate from (e.g. la-ci = MTi\^ Idt'jit,t')l3^it',t)Sait') X 
Sci{t)). In Appendices B and C we outline, in detail, the eval- 
uation of the various terms in Eq.(l33>. We will now solve 
Eq.(|33to order 0(T2). 
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E. Spin Dynamics in a Josephson Junction: An Exact Solution 

to 0{T?) 



With all of the ingredients in place, we may now solve 
EQ. (l33t to determine the spin dynamics to 0{T^). Hence- 
forth, we will examine throughout the observable "classical" 
component of the spin Sd- To make the expressions more 
appealing we will dispense with the classical "cl" subscript. 
Similarly, the action S subscript in all expectation values will 
be omitted as no time ordering subtleties appear below. We 
expand the spin as 



(34) 



Here, 5*0 is the solution to the (Larmor) problem of a sin- 
gle free spin in an external magnetic field. We computed the 
integrals borne by these zeroth order Larmor components in 
Similarly, dS{t) are the contributions borne 
These corrections 



subsection dXI C> 

by the retarded and Keldysh correlations, 
will lead to higher order contributions in (/) which are ir- 
relevant to our O(Tj^) solution. We insert Eq.( I34> into the 
equations of motion (Eqs.(I33» and retain all terms to order 
0(|Ti|2). This ti-ivially leads to 



-{SS,)-LOL{SSy) + {I^)^0, 

l{5Sy)+iOL{6S,) + {ly)=0, 



dt 



(35) 



Here, Ia=x,y.z is the a direction component of {Ici~ci + 
Iqu-ci) which was computed in the previous subsection to or- 
der 0(1 Tip). We see that the integrals / play the role of a 
driving force. Integrating, we find that 



{5 {t))^\T^\%l- cos uojt)^ 



k,p 



EkEp{Ek + EpY 
EkEpiEk + EpY 



(36) 



with 



= iTi 



k,p 



k,p 



k.p 



k,p 



; 2EkEp{Ek 


+ EpY 


. A2(^,.(0))(2^i4 




; 2EkEp{Ek 


+ Epr 




'J^j + SwLtJj) 


2EkEp{Ek - 


^EpY 


A^Symi2u;l + 





2EkEp{Ek + EpY 



(38) 



All in all, to 0{T^), the evolution of the planar spin com- 
ponents can be expressed in the format of Eas. dlOllTl with 



p{t) = COSUJLt+\Ti\ 



A2 



+ 



^ 2EkEp{Ek+Epy 
{2ijj\ + + Sw^wj) cos(ijJi + ijjj)t 

LU\ - {UJL + wj)2 

(2a;|^ + uij — SuJLUJj) cos{ujl ~ ^j)t 



(39) 



and 



/\2 

' 2EkEp{Ek + Epy 



(2a; 



- SwlWj) sin(wL - LUj)t 



(40) 



This concludes our solution for the dynamics of a spin 
in a Josephson junction. Our analysis throughout centered 
on Josephson junctions composed of s-wave superconductors 
(see our starting point EQ.( l20t ). Slightly different quanti- 
tative results appear for other pairing symmetries (allowing, 
in theory, a determination of the pairing symmetry from ob- 
servations of the spin/spin-wave dynamics and associated ef- 
fects). The deviations from simple Larmor precessions are far 
stronger for triplet (i.e. odd angular momenta) superconduc- 
tors. 



F. Physical Consequences: Josephson Nutations and Other 
Dynamical Effects 



Differentiating the equation of motion for (SS^.y) in 
Eq. (I35> and inserting the equation of motion for (SSy^x) we 
immediately obtain the equation of motion of a driven har- 
monic oscillator A simple solution yields 



{SSx{t)) = Cl cos LJLt + C2 sincjLt 

^ J COS iUnt H 2 2 ^"^) ^^"^^ 



We now discuss the physics behind our exact (to 0{T^)) 
solution. Our solution provides testimony (and to new quanti- 
tative predictions) for several, inter-related, intriguing dynam- 
ical effects. We outline these below. 

• Josephson Nutations: 

In any system harboring a continuous rotational U{1) sym- 
metry about a certain axis (z), the orbital angular momentum 
Lz is a constant of motion. Needless to say, the same triv- 
ially holds true for any spin system in which [H, Sz] = with 
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H the system Hamiltonian. In the presence of an external 
magnetic field along (or defining) the z-axis, as in the Larmor 
problem, the Hamiltonian H = —hSz commutes with Sz and 
the longitudinal magnetization (5*^ (<)) is a constant of motion. 
In our system with non local in time interactions triggered by 
superconducting pair correlations, such a conservation law no 
longer holds. Perusing Ea.(l36>. we find that the spin nutates 
above its average value. This occurrence for the S = 1/2 is 
similar to that reported in for macroscopic spin clusters 
S ^ 1. Here, however, the quantum fluctuations are pro- 
found for the 5* = 1/2 case and lead to strong deformations 
of the nutation shape. The physical engine behind the nuta- 
tions is the small time separation between the two tunneling 
electrons forming the Cooper pair As the "first" electron tun- 
nels through, it exerts a torque on the spin. A certain time later 
(of order h/A with dimensions restored) after the spin S has 
akeady revolved a small amount, the opposite spin member of 
the Cooper pair tunnels through and exerts a torque of an op- 
posite sign on the spin S. Due to the small time lag between 
the two tunneling electrons, the two opposite sign torques ex- 
erted on 5* by the two opposite sign spins of the tunneling 
singlet do not cancel and lead to a net effect. This origin is 
made evident in the retarded correlations which further 
spark a non-vanishing driving force {Id-d) along the z-axis. 
Mathematically, all of this results as the tunneling portion of 
the action contains terms which trivially do not conserve Sz ■ 
In the aftermath, this led to an effective time dependent force 
acting on Sz ■ Its form may be seen by examining the integral 
{Iz) appearing in Eq.(l35>. The latter is the z-component of 
the integrals {Id-d) and {Iqu-d) appearing in EQS.( l66ll67t . 
(Needless to say, if both members of the Cooper pair share 
the same polarization (as in triplet superconductors) then a far 
greater effect results.) 

A manifestation of the resulting dynamical effect as a con- 
sequence of these effective external forces in conventional (s- 
wave) Josephson Junctions is vividly seen in Eq. (I36> . We have 
derived similar expressions via an independent density matrix 
approach j2^. An exaggerated schematic of this effect is pro- 
vided in Fig.Q which, qualitatively, is none other than the 
standard illustration for classical rigid body nutations. We 
find that such motions now appear in the quantum arena for 
a single S* = 1/2 particle! The precise shape of our trajec- 
tories, however, is markedly different from that exhibited by 
classical rigid rotors. 

• Spin Contractions and Effective Longitudinal Fields: 
Glancing at Ea.( l67l i. the reader will see that the effective 
{Iqu-d) can be seen to dilate the spin (the uniform contribu- 
tion proportional to {S{t)) in the second equality of Eq.( I67» 
and in unison to effectively emulate a time dependent mag- 
netic field Sheff oc Cz cos <j>{t) along the z-axis in the spin 
equation of motion, d{S)/dt = ... + (S) x Sheff- Both of 
these effects were noted in j^. In EQ.( l67t . we explicitly see 
their origin. The uniform contraction is triggered by an entan- 
glement of the tunneling electrons with our S* = 1/2 particle. 
We now very briefly elaborate on the physics of this state- 
ment for the benefit of general readers. The expectation val- 
ues (S) amount to weighted sums over all possible states |?/') 




FIG. 3: The resulting spin motion on the unit sphere in the general 
case. As in the motion of classical spinning top, the spin exhibits 
undulations along the polar direction. As a consequence of entangle- 
ment with the tunneling electrons, the magnitude of the spin is not 
constant- the spin further "breaths" in and out as it nutates. 



(see Eq.( I15» . In any pure (i.e. unentangled) state of a spin- 
1/2 problem, die sum {{S oo)"^ + (Sy)"^ + {Sz^] = 1/4- the spin 
expectation values lie on the Bloch sphere. Entanglement in 
a spin- 1/2 problem such as ours is marked by a contraction, 
[{S^)^ + {Syf + {Szf] < 52 = 1/4. Any single spin expec- 
tation value within the Bloch sphere, |('(/;|S'|'0)| < S, denotes 
an expectation value computed in a multi -particle state 
which is entangled. In the case here \ip) spans the single spin 
and the tunneling electrons. Such a time dependent contrac- 
tion in the norm of (S) relative to the Bloch radius is evident 
in our exact solution of Eas.( fT0imi39ll40ll34ll36i . 

• Nonlinear planar precession: 

A notable facet of the dynamics given by the effects dis- 
cussed above are non-uniform planar precessions. We find 
that within the plane transverse to the applied field direction, 
the azimuthal angle describing the spin orientation, Lp{t) = 
ta,n~^ {{Sy{t)) I {Sx{t))) is no longer a linear in time. This ef- 
fect bears, once again, strong semblance to nutations in clas- 
sical rigid body dynamics. In the Larmor problem of a free 
spin in a magnetic field, (p{t) = lolI. In our case, the pre- 
cession about the applied field direction is no longer uniform. 
Its form is encapsulated in Eqs.( I37I3'81 or, alternatively, by 
Eas. dlOlll 1I39I40> . Once again, mathematically, the origins of 
this effect are rooted in the effective planar (xy) components 
of the effective force (/) appearing in Eq.( I35> . The explicit 
form of this effective force is given by the sum of the two in- 
tegrals evaluated in Eqs.j^ and whose origin explicitly 
lies, once again, in the same non-local in time correlations 
borne by the superconducting coiTelations. 

In summary, all of the above qualitative findings for the 
problem a single S = 1/2 spin inserted in a Josephson 
junction are made vivid in our 0{Ti) exact solution. From 
Eas. dlOll 11391 140> for the planar spin components and from 
Eqs.(l34l I36> for the longitudinal spin we clearly see how all 
of these effects come into play. 
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VI. RETARDED CORRELATIONS IN GENERAL SPIN S 
DYNAMICS IN A JOSEPHSON JUNCTION 

The equation of motion, Eq.(|9}, is valid for all spins 5'. 
Much of our formalism follows with no change. We now ex- 
amine the integrals {!) in the general spin S problem. 

We find that for a spin of size S, the integral {Iqu-d) un- 
dergoes no change relative to its S = 1/2 form- Ea.( l67t re- 
mains the same. The associated physics fleshed out in the 
second equality of Ea.( l67t which was described in the previ- 
ous section (spin contractions and the presence of an effective 
longitudinal field) undergoes no change for the general spin S 
case. 

Next, we evaluate {Id-ci)- For large spins, 5* ^ 1, the 
product {S{t) X S{t')) is well approximated by the vector 
product of averages {S{t)) x {S{t')). Then, the approach of 
Ref. 13] is well justified and we can obtain the Josephson nu- 
tations of a big spin. For any S and t' > twe obtain 

{Sd{t') X Sd{t))s 
= l{[{Sy{t), S.m+ [cosuLif -t)-l] 

- {S4t),S,it)}+ sincuLit' -t)]e^ 

+ [{SAt), S,{t)} +[1- cos LUL{t' -t)] 

- {Sy{t),S,it)}+ sinLULit' - t)] iy 

+ {2Sl{t) + 2Sl{t))smujL{t' -t)e,), (41) 

where {...}+ denotes an anticommutator In the S = 1/2 
case, the integral stemming from these vectorial product cor- 
relations is parallel to the z-axis, {Ici~ci) \ \&z- For spins of size 
5 > 1/2, however, as we see from Ea. (l41> . the planar (a;, y) 
components also come to the fore and lead to retarding corre- 
lation (/3^) effects in Ea.(l33>. Furthermore, the magnitude of 
the driving force {Ici~ci;z) along the z-axis, much unlike the 
= 1/2 case is time dependent. 

For general spins of size S > 1/2, both retarded (/9^) and 
Keldysh {(3^^) are non-zero along any spin direction. All of 
the effects discussed in subsection dV F> are present. 

It is noteworthy to discuss the scaling of all terms with the 
spin size S. As evident from Ea.(l41>. the integral {Id-d) 
spawned by retarded correlations scales as S^. Similarly, as 
seen from Ea. (l67> . whose form holds for arbitrary S, the ef- 
fective driving force {Iqu-d) generated by Keldysh correla- 
tions {(3^) scales as S, i.e. {Iqu-d) oi S. Thus, for large 
spins S ^ 1, the retarded contributions overwhelm stochastic 
Keldysh contributions. In the classical limit, S oo, only 
the retarded contributions remain iHl . 



VII. SPIN TRIGGERED AC EFFECTS 

Thus far our discussion centered on a Josephson junction 
for a time independent potential difference V (dc voltage bias) 
between the two superconducting leads for which (f){t) — ujjt 
with u!j = 2eV. 



We now briefly sketch matters for an ac voltage bias 
wherein the potential drop is oscillatory in time and the cor- 
responding phase difference is (p{t) = Asin^lt. To make the 
physics more transparent, we omit any dc contributions to the 
voltage (and thus linear in time contributions to the phase). 
This serves as a caricature of rf driven Josephson junctions 
known to exhibit the famous Shapiro steps lIZTIl . 

The setup is given by Fig. |2] for a spin S = 1/2 particle 
yet now with an ac voltage applied across the junction. In the 
sections that follow, we will resume our central focus on the 
constant voltage drop case, ^(t) = cujt. Only in this short 
section do we analyze an applied ac voltage bias. 

The calculations for the ac voltage bias case parallel the 
analysis of the previous sections. First, we express all terms 
by pure harmonics. This is readily achieved by relying on the 
identity 



(42) 



with { J„(C)} Bessel functions. The factor t') of Ea.( l32l i 
and thereafter now becomes 

i(t,i') =^^(^)^m(^)cos[l](nt + m<')]. (43) 



The analog of Ea. (l66> for the ac voltage bias case is 

{Id~d)s - 

2mD,ujL sin r2(n + m)t 



(44) 



Further resonant (delta function) terms make an appearance 
for m 1. 

Similarly, the analog of Ea. (l67> reads 

n.rn,k,p 

X [{2mQ{Sx{t)) sin(n + m)flt — {Sy{t))LJL cos(n + m)Qt)ex 
+ {2mi}{Sy{t)) sm{n + m)Qt + {Sx{t))uJL cos(n + ■m)Qt)ey 

+2'mn{Sz{0)) sin(7i + m)rite^]. 

To OdTi p), the nutations are given by 



2mCl , lul 
X— ^^i- 



{Ek + Epy'Ek+Ep 



1 — COS fl(n + m)t 
X ^ L, (45) 

n + m 

Higher order effects further enhance this response. Ea.(l45> is 
the ac voltage bias analog of Eq.(l36> for the dc voltage bias 
case. The seminal feature of our results is the existence of fre- 
quencies in the spin dynamics of all integer multiples of the 
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voltage bias driving ac frequency 51. As the spin alters (via 
back-action effects) the tunneling supercurrent, the supercur- 
rent will exhibit oscillations at all frequencies ujr = r^l with 
r an integer Extending the results of |3] to this problem, the 
supercurrent 



(1(1)) =sin0(t) 2Tr^ep^A{\To 



\Ti\') 



+4e\Ti\'p'h(S,{t)) , 



(46) 



where p is the spin density of states within the leads, with 
the spin given by Eqs.(I34| 145 > with the Larmor {Sz{t))o = 

(^.(0)>. 




FIG. 4: A ferromagnetic slab inserted between two superconducting 
leads. The entire system is subjected to a weak external magnetic 
field B. A schematic of the precessing spins is shown. 



VIII. FERROMAGNETS IN JOSEPHSON JUNCTIONS 

We now investigate what transpires when ferromagnets (in- 
stead of a single spin) are immersed between two s-wave su- 
perconductors with a dc bias voltage applied across the junc- 
tion (as illustrated in Fig.Q). As in the single spin problem, 
the full problem involves both the back-action of the spin on 
the phase of the superconductors (ignored here) and the spin 
dynamics sparked by the tunneling current (which we focus 
on below). Further, for extended junctions, phasons naturally 
appear. In what follows, we assume that the phases of the two 
superconductors surrounding a single magnetic slab have a 
spatially uniform phase difference (/)(t). The tunneling action 
amounts to a sum over individual tunneling actions through 
each of the individual spins labeled by their sites r, 

Stunnel — 



+ JdtJ dt'Ttl3'^{t,t')S,u{r,t) ■ S,u{?,t')j{t,t'). 

For ferromagnetic spin chains/planes with arbitrary exchange 
constants J{r,r'), and scaled external magnetic field h, the 
exact equation of motion reads 

= (^^^^^ +hx Scl{r) + Icl-chr + Iqu-cl-r 



dt 



- ^ J(r, 7^)5,z(r', t) X 5,z(r, t))s. (47) 



It is hard not to notice a resemblance between the single 
spin problem (Eq.j33t) and the problem of the ferromagnet 
(Eq.(l47h. Indeed, as we will shortly demonstrate the spin 
wave dynamics in the ferromagnet within a Josephson junc- 
tion bears much in common with the single spin problem with 
the proviso that the various ferromagnetic spin waves feel 
an effective momentum dependent magnetic field of strength 
h^ff ^ h + S[J{k) — J(0)] with J{k) the Fourier transform 
of the two spin interaction J(r, f*). 

The solution proceeds much the same as for the single spin 
problem. Henceforth, we discuss the qualitative physics ex- 
pected. Unlike the precise solutions presented till now, what 



follows is a quick qualitative sketch by wayof an analogy. An 
exact solution will be detailed elsewhere y31 . 

Transforming from spin variables to bosonic operators 
(6(r)) at all lattice sites r, IH 

= b\f)V2S, 



(48) 



S-if) = [6(r) - —b^r}bir}bir}W2S, 
= -S + b^{r)b{r). 



Sans the 0(1 Tip) tunneling part of the action, the action is 
quadratic in the bosonic operators and is readily diagonalized 
in q space. We find that the free part of the action 



Sq^- dt 



d'^q 



{2n) 



^{S[J{q) ~ J(0)] + h} 



Xb*iq)b{q), (49) 



with d the dimension of the inserted magnet. (As the problem 
is ferromagnetic, J(0) = min^ J((j)}). Comparing this ac- 
tion to the one appearing in the single spin problem, we find 
that to Gaussian order the spin-wave problem is identical to 
the dynamics of a single spin with the replacement 



h^Kff{q)^{S[Jiq)^jm + h}. 



(50) 



The quadratic contribution of the 0(|rip) portion of the ac- 
tion involving non-local in time correlations has precisely the 
same form for both the single spin problem and for each mode 
q of the spin- wave problem. Thus, the quadratic in fe, OdTi p) 
corrections to the spin dynamics are given by Eqs. j36l37l with 
the replacement of Eq.dSOIl. 

For instance, the above analogy suggests that that the net 
ferromagnetic moment variation in 5 = 1/2 ferromagnets is 



SM 
~V~ 



|Ti|2(l-cosa;,7t)[^ 



k.p 



EkEpiEk + EpY 



A2 



M 



V ^EkEpiEk + EpY' 



(51) 



with V the volume of the magnet and M its magnetization 
sans the supercurrent. Alternatively, the analysis may parallel 
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the derivation of the previous sections word for word while 
taking the unperturbed solution (the analogue of the Larmor 
solution of Eq. (I61» to be a spin wave and computing all cor- 
rections to ©(Tf). 

In the continuum limit. 



mo 



(52) 



with mo = S/ V (where v is the volume per site), the magne- 
tization density of the ground state, and ps the spin stiffness. 

This h — > hef f (q) correspondence applies to any property 
inherited by the single spin dynamics in the Josephson junc- 
tion. In particular, in Q it was beautifully shown how spin 
dynamics may alter the super-current in the Junction. The 
current may be computed by the likes of Eg. (127k Extending 
these results to a ferromagnet inserted in a Josephson junction 
by the correspondence of Eq. (I50> . we find that the new spin 
wave dynamics leads to the supercurrent. 



{I {t)) = sin (f>{t) 



^[2.VA(|ToP-J|r,n 

+4e|Ti|V/ie//(g1(^.(g1)l,(53) 



with p the spin density of states within the leads. A match- 
ing of the Josephson and spin frequencies (such as present 
here for variations in the low temperature magnetic moment 
(see Eas.( l36l5'll )) leads to a DC signal; additional harmon- 
ics further appear. We emphasize that in the above we com- 
pared only the Gaussian portion in the Bose fields. Higher 
order (non-Gaussian) terms originating from Eq.( I48> as well 
as phasons alter the natural correspondence of Eq.(I50>. A full 
discussion of these issues will be detailed elsewhere 12^. 



integrals to low orders was replaced by an "adiabatic" ap- 
proximation wherein the slow dynamics of the spin vis a 
vis electronic processes was explicitly incorporated, S{t') ~ 
S{t) + {t' -t){dS/dt). 

(ii) The ("classical") large S limit allowed us to omit many 
instances of Squ in the equations of motion and only (3^ re- 
lated contributions in the tunneling action were consequential. 
Furthermore, as briefly alluded to earlier, in this limit, the av- 
erage of the vectorial product {S{t') x S{t)) is equal to the 
product of the averages {S{t')) x {S{t)). Correspondingly, 
any expectation value braces may be omitted. Thus, we may 
replace any expectation value {A) by A itself. 

The advantage of this method is that furnishes an elegant 
non-perturbative closed form solution for the spin dynam- 
ics. We will not repeat the results for the single spin cluster 
(5 S> 1) problem here and rather refer the reader to jsl. We 
now briefly comment on applications of this method to other 
systems. 

To OdTi p), the spin wave dynamics in ferromagnets may 
be attained via the substitution of Ea.( l50> . Equivalently, the 
spin wave equations of motion may also be determined di- 
rectly when applying the adiabatic approximation on Eq.(l47r 
We then find 



dt 



h X Sciin) + ^ JijSci[rj,t) X Sci{ri,t) 



dS 

-kSci X — — smLUjt = (54) 
dt 



with ^ ^ Ek.p^^liEk +E,- eVr' {E, 
eV)^'^]. The appropriate spin wave equation is 



k + Ep 



IX. OTHER GEOMETRIES 

If phason contributions are neglected, then by a trivial 
change of geometry all of our results thus far, will apply 
for other systems as well. For instance, by replacing one of 
the superconducting leads by a surface, the resulting system 
may emulate a superconducting tip coupled to superconduct- 
ing surface through a single spin or a ferromagnet. Here, all 
of the results of Sections dVIVllU for the spin dynamics and 
tunneling current hold. 

Similarly, by replacing both superconducting leads by sur- 
faces and examining a magnetic layer inserted in between, 
the resultant system looks much alike a layered supercon- 
ducting/magnetic system. In this system, the results of 
Section fVllH apply. 



X. LARGE S ADIABATIC APPROXIMATIONS 

Thus far we studied the dynamics of single spins and of fer- 
romagnets. In 0], the large S limit of the single spin problem 
was studied. In that work, several approximations were made: 

(i) The perturbative approach that we employed in the cur- 
rent article which allows an exact evaluation of all pertinent 



dm 

dt 



z[h + S{Jiq)~JiO)}]b{q) 



+Kdtb{q, t) sinujjt. 



(55) 



The solution to Eq.(l55> is 



= 6(q',0)cxp[ 



2i{S{J{q)-J{0)} + h) 



x{tan-i(^==)-tan-i( ,--4 1], (56) 



which is quite different from the standard spin-wave evolution 
in a magnet outside a Josephson junction. The key feature 
is a nonuniform evolution of each spin-wave. Similar to the 
azimuthal precession of a single spin, the planar components 
Sx,y precess as the real and imaginary parts of exp[i(^(t)] with 
a nonlinear (^(t). Thermodynamic quantities computed via the 
corrected bosonic spin- wave dispersion exhibit corrections. 

Similarly, we may examine the adiabatic large S equations 
of motion for an antiferromagnetic spin chain oriented along 
the z-axis in a Josephson junction (just as in Fig.Q yet now 
with a single antiferromagnetic spin chain replacing the ferro- 
magnetic slab in an = background). We then find that that 

the staggered spin. Si = (— l)*S'i (with the integer i the spin 
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site location along the chain) satisfies 

+KdtSci (t) sin ujjt, (57) 

where □ = ^ [9^ — -^d^]. Here, g = 2/S and the spin wave 
velocity Vg ~ 2a J S, with a is the lattice constant. The role of 
the supercurrent as an effective driving term is evident in the 
last line of Eq.j57>. 

XI. CONCLUSIONS 

In conclusion, our work addresses new dynamical effects 
exhibited by spins in Josephson Junctions. En route, many 
features (general and specific) were found: 

(1) We derived the exact equation of motion for spin sys- 
tems on Keldysh contours. 

(2) The S = 1/2 spin dynamics of a single spin in a Joseph- 
son junction was investigated and a perturbative solution was 
given. Spin-1/2 Josephson nutations are predicted. 

(3) Spin dynamically triggered ac effects are predicted. 

(4) The spin wave dynamics of a ferromagnet in between 
two superconducting leads was investigated. We predict non- 
trivial spin wave dynamics as well as new manifestations of 
this dynamics (most notably in the supercurrent). 

(5) Large S expressions were discussed for ferromagnetic 
slabs and antiferromagnetic spin chains in a Josephson junc- 
tion. 
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A. Appendix A: Detection 

The non-trivial spin-wave and associated supercurrent in 
Josephson junctions containing ferromagnets (section IVllH 
may be seen more readily seen than those of single spins. The 
spin dynamics may be discerned by measuring the magneti- 
zation of the ferromagnetic slab as a function of time (as sug- 
gested by Eq.(I51» as well as by monitoring the supercurrent 
(given by Eq.( I53» . Other techniques may involve standard 
measurements of microwave radiation from the junction (and 
backaction effects). The magnitudes of these effects will be 
studied elsewhere l29il . 

We now briefly review a detection scheme discussed in 
01,113] for the Josephson nutations for the 5* ^ 1 limit of 
the general spin S results of Section fvU . This corresponds to 
a single magnetic cluster 

As it moves, the spin cluster magnetic moment generates a 
time-dependent magnetic field, 6B{r, t) = ^ [3f(r • Tn{t)) — 
r'^rh{t)]/r^. This small field is superimposed against the con- 
stant external field background (B). In the above, r is the 




FIG. 5: A SQUID-based detection scheme. The SQUID monitors 
the magnetic field produced by the magnetic cluster in one of the 
junctions. 



position relative to the spin and m(t) is the magnetic moment 
of the spin. A ferromagnetic cluster of spin S = 100 gen- 
erates a detectable magnetic field SB ^ 10^^° T at a dis- 
tance of a micron away from the spin. For a SQUID loop of 
micron dimensions located at that position, the ensuing flux 
variation (5$ ^ 10~^$o (with $o = hc/e the flux quan- 
tum) are within reach of modern SQUIDs. In such a setup, 
with Ti/Tq ^ 0.1, the typical critical Josephson current is 

4°^ - 10 pA, |A| = 1 meV, and eV - lO^^jAI, we find 
that the relative corrections SS/S ^ 0.1. The spin compo- 
nents orthogonal to B vary, to 0{Ti ), with Fourier compo- 
nents at frequencies jw^ ± loj] {ujl = SPbB), leading to a 
observable signal in the magnetic field B + SB. For a field 
B ~ 200 G,ujl ^ 560 MHz, and a new side band will appear 
at \ujl — ujj\, whose magnitude may be tuned to 10-1 00 MHz . 
This measurable frequency is easily distinguished from the 
Larmor frequency ui l ■ 

The efficiency of this detection scheme may be enhanced by 
embedding the spin in one of the Josephson junction arms of 
the SQUID itself. Such a setup is illustrated in Fig. |5] The 
Josephson junction harboring the spin is employed in both 
triggering the nutations and, along with the second junction 
of the SQUID, in the detection of the resulting nutations. 



B. Appendix B: Time Ordering Along the Keldysh Contour 

In averaging within the path integral formalism, we imme- 
diately attain time ordered averages. In interchanging the or- 
der of the spins (if necessary) in the vectorial product upon 
time ordering within the path integral CPi formulation a 
change of sign is incurred |31]. We now go over, in some de- 
tail, time ordering within the Keldysh framework. As the time 
ordering is performed along the Keldysh contour, we will de- 
note it by Tk as we have done in deriving the effective action 
of subsection dV C> . 

Consider the third term in Eq.( I33> . Upon time ordering, we 
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find that 

(V(t') X Sci{t))s = {TK[S,u{t') X S,im 

= {T^Kll^i^upit') — Sdown{t')) 

^{Sup{t)+Sdo^n{t))]). (58) 

Due to the form of the Keldysh contour (see Fig.Q), irrespec- 
tive of the values of t and t' , Sup{t) always appears before 
Sdown{t'). Similarly, fort > t', Sup{t) appears after Sup{t') 
while Sdown{t) appears before Sdown{t')- With this informa- 
tion at hand, Ea.(l58> leads to 

{TK[Sgu{t') X S,im 
= -e{t - t')[{S{t') X S{t)) + {Sit) X Sit'))]. (59) 

The expectation values on the right are the usual operator ex- 
pectation values. Here, we disposed of the up/down indices 
once we took care of time ordering. The up/down labels 
merely serve as mnemonics for this time ordering along the 
Keldysh contour 

Similarly, we find that 

{TK[Scl{t') X Sclit)]) = 



-[eit'-t)i{sit')xsit)~sit)xsit'))) 

+0it - t')i{Sit) X S{t') - Sit') X Sit)))]. 



(60) 



By the same token, (T^ [5g„(t) x Squit')]) = 0. 

As will become clear shortly, in the solution of Eas. (l33> to 
order ©(Tf ), we will need the spin-spin expectation values 
of the usual Larmor problem (i.e. a single spin in a magnetic 
field sans any supercurrent). Here, 

Sxit) — SxiO) coaojLt + SyiO) sinojLt 
Syit) — SyiO) cosu) i^t ~ S xiO) sin uj^t, 

S^t) = '^.(O), (61) 

with the external magnetic field oriented along the positive z 
axis and lul = |/^| the Larmor frequency. 

Next, we invoke this solution to compute the various expec- 
tation values within the Larmor problem (i.e. to order ©(r}*)). 
We find that 

{Sguit') X Sclit))s 

= -9it - t')[i Im{{Sit') X Sit))} 

= -lOit - t'){[(5,(t))(l + COSUJLit' - t)) 

+ {Syit))s\TLUJLit' ~t)]e^ 

+ [{Syit))il + COSLJLit' ~t)) 

-{S^it))smujLit' 't)]ey 
+2{S,it))cosLJLit' -t)e,}. (62) 

Similarly, for the 5 = 1/2 problem, 

{Sciit') X Sciit))s 
= Re{{Sit') X Sit))} 



In EQS. J62l6'3t . we vividly see that upon time ordering 
along the Keldysh contour, the non-vanishing spin cross prod- 
ucts become simply related to the imaginary and real parts of 

{Sit') X Sit)). 



C. Appendic C: Evaluation of Integrals 

We are now ready for the evaluation of the various integrals 
/ that appear in Eq.(|33} to order 0(rf ). 

We start with Ici~ci- Inserting Eq.(l63ll and Eq.(l30ll into 
Eq.(l33> we find upon invoking the relation (pit) ~ ujt 

{Icl-cl)s = 

4|Ti|2e, [ dt'jit,t')/3''it,t'){S,iit') X S,iit))s 
= ^]T,]'e^J2l^ I dt'eit-t')cosiUj-^^' 

X siiiULit - t') sm[iEk + Ep)it - t')]. (64) 



Before evaluating Eq.(I64> exactly, we illustrate what an- 
swer is anticipated. The underlying observation of this adi- 
abatic approach is that, as a consequence of ojl.j <C Ek^p 
the spin dynamics is far slower than that of electronic pro- 
cesses. Thus, in integrals involving both spin and electronic 
degrees of freedom, we may regard the spin as nearly sta- 
tionary and approximate Sit') ~ Sit) + {t' — t)idS/dt). 
This physically transparent approximation was invoked in 
0]. Employing this approximation here we anticipate that 
{Ici-ci)s — CiCzSinwii where Ci is, up to trivial prefac- 
tors, given by dx [x^ (3^'ix)]. Such an anticipation is not 
far off the mark. 

Next, we exactly evaluate Eq. J64> by rewriting products of 
trigonometric functions as sums and consequently employing 
the identities 



dx cos ax = 7r(5(a), 



dx sin ax = 7r(5(a) 



(65) 



In the integrals of interest, x assumes the role of {t' — t). As 
the applied magnetic field and voltage are far lower than elec- 
tronic energy scales, lul.j ^ Ek.p, we find that all resonances 
signaled by the delta functions are physically unaccessible and 
our expressions undergo further simplifications. Retaining the 
leading order terms in OiujL,j /{Ek + Ep)) we arrive at 



{I 



E 



EkEpiEk + Ep)- 



■sinujjt. (66) 



■ smtOLit' — t)ez 



(63) 



Thus, the form anticipated by the adiabatic approximation is 
correct if Ci = -|TiPE.,,^;:^feffcF- 
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Similarly, by inserting Eos. ( 13 11621 1 into Eq.( I33> and invok- 
ing Eqs.( I65> . we find 

{I,u-ci)s = -iTil'Yl^^^^^yi^TE;^ 

x[i{Sxit))ujj amcjjt - {Sy{t))bJL cosujjt)ex 
+{{^yi't))'^J sinujjt + {Sx{t))uJL cosujjt)ey 
+ {Sz{0))ujj ainujjtez] 
= C2lLUj{S{t)} sinc^it) + LJL{e, x (S {t))) cos c^{t)] , (67) 



with the constant Ca = -\Ti\^Ek.p ig.gjgl+g,)^ ^nd 
<j>{t) ~ LUjt the superconducting phase difference across the 
junction. The last line of Ea A61l has a very physically sug- 
gestive meaning regarding spin contractions and an effective 
longitudinal magnetic field- items which we will expand on in 
Section (IVF> . Our expressions (Eqs. (I66I67» above are exact 
to lowest order in Ti and the ratios {ujl.j / Ek,p)- 

Finally, the integral {Iqu-qu) = identically by virtue of a 
vanishing {Squ{t') x Squ{t))s = 0. 
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